Electronic shell effects and the stability of alkali nanowires 
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Experimental conductance histograms for Na nanowires are analyzed in detail and compared 
to recent theoretical results on the stability of cylindrical and elliptical nanowires, using the free- 
electron model. We find a one-to-one correspondence between the peaks in the histograms and the 
most stable nanowire geometries, indicating that several of the commonly observed nanowires have 
elliptical cross sections. 
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INTRODUCTION 

In the past eight years, experimental research on metal 
nanowires has burgeoned, together with theoretical work 
describing these nanoscale systems in which quantum ef- 
fects play an important role (see Ref. [1] for a recent 
review). Metal nanowires exhibit a number of interest- 
ing properties: their electrical conductance is carried by 
a well-defined number of quantum channels [2, 3, 4, 5], 
their shot-noise is suppressed by the Pauli principle [6], 
and they are remarkably strong and stable [7] . The atten- 
tion to this research field is accentuated by the ongoing 
miniaturization in the microelectronic industry and the 
need to understand what modified properties may be met 
upon further size reduction. 

Crucial for answering the question of future appli- 
cations is knowledge about the stability of nanowires. 
An important contribution to the energetics of metal 
nanowires comes from the subband structure of the con- 
duction electrons [8, 9, 10]. Indeed, a linear stability 
analysis [11, 12, 13, 14] shows that electronic shell effects 
can stabilize and therefore favor certain geometries. The 
effects of shell filling on the abundance spectra of metal 
clusters have been known for years. They explain the 
existence of clusters with "magic numbers," correspond- 
ing to full electronic shells, which are observed more fre- 
quently than others [15]. More recently, electronic shell 
and supershell structures have been reported for alkali 
metal nanowires [16, 17, 18]. Theoretically, these effects 
can be understood within a stability analysis including 
wires with cylindrical and elliptical cross sections [14]. 

This paper is organized as follows: An overview on the 
theoretical framework is presented in the next section, 
including the discussion of the nanoscale free-electron 
model. Then experiments measuring conductance his- 
tograms of alkali metal nanowires are presented, followed 
by a more detailed analysis of the experimental data. It 
is shown that theory allows prediction of the shapes of 



wires produced in experiments. A short summary is given 
in the last section. 

THEORY 

One possibility to describe nanoscale contacts is via 
molecular dynamics simulations [19, 20, 21] which use 
short-ranged interatomic potentials suitable to describe 
the bulk properties of metals. However, such an approach 
appears problematic when applied to metal nanowires, 
in which electron-shell effects [16, 17] due to transverse 
confinement are likely to be important. An alternative 
approach, motivated by the success of the jellium approx- 
imation in explaining the energetics of ultrasmall metal 
clusters [15, 22], is a nanoscale generalization of the free- 
electron model, developed by Stafford et al. [8]. This 
approach allows for analytical calculations and was suc- 
cessfully used in recent years to understand many prop- 
erties of metal nanowires. 

Nanoscale Free-Electron Model 

The free-electron model of a nanowire [8] treats the 
electrons as a non-interacting Fermi gas kept within the 
wire by a confining potential, which is approximated in 
terms of hard-wall boundary conditions. The discrete 
ions are replaced by a constant positive jellium back- 
ground charge. This model requires good charge screen- 
ing in the metal and an almost spherical Fermi-surface. 
It is especially suitable for alkali metals, but can also be 
applied to other monovalent metals such as gold. 

A nanowire connecting two macroscopic electrodes is 
an open quantum system, for which the Schrodinger 
equation is most naturally formulated as a scattering 
problem. The fundamental quantity describing the prop- 
erties of the system is the scattering matrix S(E) con- 
necting incoming and outgoing asymptotic states of con- 
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duction electrons in the electrodes. Transport properties 
can be expressed in terms of the submatrix S12 describing 
transmission through the wire. For example, the electri- 
cal conductance G is given by the Landauer formula [23] 
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where f(E) is the Fermi-Dirac distribution function. 
Similarly, thermodynamic properties can be expressed in 
terms of the scattering matrix through the electronic den- 
sity of states (DOS) D(E) 
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from which the relevant thermodynamic potential for an 
open system, namely the grand canonical potential £1, is 
obtained as 



fl = -k B T J dE D(E) In 
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where ks is the Boltzmann constant, T is the temper- 
ature and 11 is the chemical potential specified by the 
macroscopic electrodes. Eqs. (1) through (3) include a 
factor of 2 for spin degeneracy. 

Many fundamental phenomena of metal nanowires 
have been understood in terms of this nanoscale free- 
electron model. In particular, theoretical results on con- 
ductance quantization [24, 25, 26], cohesion [8, 26] and 
current noise [27] are in good agreement with experi- 
ments. At the mean-field level, electron-electron inter- 
actions can also be included in the model in a straight- 
forward way [26, 28, 29]. 

In this approach, an atomic-scale contact between two 
macroscopic electrodes can be considered as a waveguide 
for conduction electrons (which are responsible for both 
electrical conduction and cohesion in simple metals) [8]. 
Each quantized mode transmitted through the contact 
contributes one quantum Go = 2e 2 /h to its conductance, 
and acts as a chemical bond, dclocalized through the con- 
tact, thus contributing to the cohesion of the wire. 



Stability and Shell Structure 

Metal nanowires are found to be remarkably stable, 
although one might expect them to break up because of 
surface tension. Such an instability arises from classical 
continuum mechanics and is known as the Raylcigh in- 
stability [30]. The key to the stability of nanowires is 
quantum corrections to the classical stability analysis. 

For wires aligned along the z-axis, the geometry is 
given by a function R(z, (p) when using cylindrical coordi- 
nates. The linear stability of a nanowire with a given ge- 
ometry can be analyzed by calculating the energy change 
<50 caused by a small deformation 5R(z, ip) of the wire 



geometry. A stable realization of a nanowire requires an 
increase of fi for all possible small deformations. A con- 
straint on the deformation arises from the fact that, de- 
pending on material parameters, the deformed wire tries 
to find a compromise between a volume conserving defor- 
mation and one ensuring electroneutrality [28] , leading to 
the constraint 
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where V is the volume of the wire, S its surface area, and 
kp = 2-k/Xp the Fermi wavevector of the electrons. The 
parameter 77 can be adjusted so as to fix the value of the 
effective surface tension to the material-specific value. In 
particular, r\ — corresponds to a constant-volume con- 
straint, and r\ = 1 is the constraint of constant Weyl 
charge. In this paper, the value 77 = 0.93 is used, corre- 
sponding to the surface tension of Na (0.22 N/m [31].) 

The stability was first analyzed for cylindrical geome- 
tries and axisymmctric deformations, since these arc the 
only destabilizing deformations in the classical calcula- 
tion. In a scmiclassical approach [11, 12], the grand 
canonical potential is divided into a smooth term de- 
pending on classical geometric quantities (volume, sur- 
face, mean curvature), referred to as the Weyl con- 
tribution [32], and a fluctuating quantum term. This 
quantum potential is calculated using a Gutzwiller- 
type trace formula [32, 33]. One finds that for cer- 
tain radii electronic shell effects suppress the classical 
Rayleigh instability at low temperatures. This leads to 
intervals of stable ("magic") radii, resulting in a series 
of stable cylindrical wires with conductance G/Gq e 
{1, 3, 5, 6, 8, 10, 12, 14, 17, 23, . . .}. Although all of those 
wires are predicted to be metastable, some of them 
(G/Go = 5, 14, . . .) are expected to be barely stable, and 
thus unlikely to be observable. Since the semiclassical 
theory depends only trivially on the length of the wire, a 
puzzling conclusion is that nanowires with a magic radius 
remain stable whatever their length. 

A full quantum calculation [13] shows that the suppres- 
sion of the Rayleigh instability by electron shell effects is 
supplemented by an interplay between the Rayleigh and a 
novel Peierls-type instability missing in the semiclassical 
approximation. In fact, this latter quantum mechanical 
instability limits the maximal length of stable nanowires. 
For lengths of order 10-1000 nm, the quantum calcula- 
tion mainly confirms the set of magic conductance val- 
ues mentioned above, but their corresponding intervals 
of stable radius values are reduced. Wires with a con- 
ductance of 5Go and 14Go are unstable and therefore 
missing. 

Axial symmetry implies characteristic gaps in the se- 
quence of stable nanowires, which is not fully consistent 
with the experimentally observed nearly perfect period- 
icity of the peak positions as a function of \fG. The 
deviations can be accounted for neither by the inclusion 
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FIG. 1: Calculated energetically stable cylindrical and ellip- 
tical nanowires as a function of temperature (in units of the 
Fermi temperature Tp). The surface tension is adjusted to a 
value of 0.22 N/m (77 = 0.93), corresponding to Na [31]. The 
aspect ratio e is coded via the gray scale, and hatched areas 
indicate multistability 



EXPERIMENTS 

Using remarkably simple experimental techniques, it 
is possible to gently break a metallic contact and thus 
form conducting nanowires. Besides the use of scanning 
tunneling microscopes (STM), the mechanically control- 
lable break junction (MCBJ) is the most commonly used 
tool for the study of mechanical and transport prop- 
erties of metal nanowires and point contacts. By in- 
denting one electrode into another and then separating 
them, a stepwise decrease in electrical conductance is ob- 
served, down to the breakpoint when reaching a size of 
a single atom in cross section. In addition to the con- 
ductance, the STM and MCBJ techniques were used to 
study a variety of properties of metal nanowires, i.e. 
cohesive force [37], conductance fluctuations [38], ther- 
mopower [39], and shot noise [6]. The highly reactive 
alkali metals have been studied using a slightly modified 
MCBJ technique [4]. 



of disorder [34], nor by the use of more elaborate self- 
consistent jellium models [35, 36]. Since gaps in the se- 
quence of cylindrical nanowires arise from a degeneracy 
of conductance channels, it is natural to assume that a 
Jahn- Teller deformation breaking the symmetry can lead 
to more stable deformed configurations. 

In fact, the stability of nanowires with elliptic cross- 
sections was recently examined theoretically [14] and sta- 
ble elliptical nanowires were found. This extended anal- 
ysis confirms the stability of the cylindrical wires found 
in the quantum calculation, even with respect to Jahn- 
Teller distortions. This includes the fact that cylindrical 
wires with a conductance of 5Go and 14Gn are energet- 
ically unstable. Secondly, a number of stable elliptical 
nanowires with an aspect ratio e (which is defined as the 
ratio of the two major semi-axes of the ellipse) greater 
than one are found. The most stable elliptic wires have 
conductance values of G/Go € {2, 5, 9, ...}. 

Fig. 1 shows the corresponding stability diagram in 
the low conductance range as a function of temperature. 
Shaded areas indicate stability and the gray scale repre- 
sents the value of the aspect ratio e. The x-axis is given 
by the corrected Sharvin conductance Gs of the wire [1] . 
The elliptical wire with conductance 9 Go is expected to 
be much more stable at low temperature than the neigh- 
boring cylindrical wires at 8 Go and 10 Go, because it is 
stable up to a much higher temperature. 

Extending this analysis to higher conductance values 
and extracting the most stable configurations [14] leads 
to a sequence of stable cylindrical and elliptical nanowires 
that allows for a satisfying interpretation of the experi- 
ments on shell and supershell structure [16, 17, 18]. 



Conductance histograms 

Nanowires fabricated by the methods discussed above 
show a wide variety of behaviours, with the atomic struc- 
ture of the contacts playing an important role. For ex- 
ample, each scan of the dependence of conductance G 
on the elongation d is individual in detail as the atomic 
configuration of each contact may be widely different. 
Generally, the experiment doesn't allow for a detailed 
knowledge of the atomic-scale structure of the contacts. 
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FIG. 2: Experimental conductance histograms for Na at dif- 
ferent temperatures. The curves are vertically offset for clar- 
ity. The black arrows mark the approximate positions of the 
peaks that we consider in our analysis. (Data taken from 
Ref. [16].) 



Nevertheless it is possible to extract information about 
possible quantization effects by performing a statistical 
analysis. Many individual conductance versus displace- 
ment curves together produce a histogram of the proba- 
bility for observing a given conductance value, which is 
quite reproducible for a given metal and for fixed exper- 
imental parameters. 

An important subject for study is the alkali metals, 
since they are nearly-free electron metals and therefore 
most closely approach the predictions of the nanoscale 
free-electron model, described in the previous section. 
Fig. 2 shows conductance histograms for sodium recorded 
at different temperatures (taken from Ref. [16]). The 
curves are vertically offset for clarity and the ampli- 
tude has been normalized by the area under each graph. 
Each histogram is constructed from 1000-2000 individual 
scans. A number of peaks with very different intensity 
can be identified: Sharp peaks can be recognized at con- 
ductance values near 1, 3, 5 and 6 Go, while rather wide 
maxima are found near 9, 12 and 17 Go- In addition, 
there is a small peak near 2 Go and a very weak struc- 
ture near 4 Go- Results similar to those for Na have been 
obtained for Li and K. 



FURTHER ANALYSIS OF NA HISTOGRAMS 

In order to obtain more information about the stable 
geometries that are realized in experiments, we have an- 
alyzed in detail the conductance histograms for sodium 
recorded in the temperature range of 20-70K (Fig. 2). 
We fitted each data set with a function consisting of a 
sum of peaks plus a background. The conductance peaks 
were taken to be of lorentzian shape with height, width 
and peak position as fitting parameters. We have cho- 
sen to take nine peaks, corresponding to those that are 
clearly distinguishable in the data set. Their approxi- 
mate positions are marked with black arrows in Fig. 2 
and they are labeled with capital letters. This choice 
is not unique, as the fit can always be improved by in- 
cluding additional broad peaks of lower intensities. How- 
ever, the parameters of additional peaks can not be de- 
termined with high confidence, especially for noisy curves 
at higher temperatures. Fig. 3 shows the fit for the data 
at T = 30K and the inset shows a fit of the region near 
9Go, including two additional peaks of lower intensity, 
for comparison. 

The values of the fitting parameters now allow to ana- 
lyze the development and the contributions of the differ- 
ent conductance peaks as a function of temperature. The 
peak positions stay almost constant with maximum fluc- 
tuations of three percent within the temperature range 
from 20-70K. All the peaks have a finite (temperature 
dependent) width, and their positions arc shifted down- 
ward compared to the theoretical quantized conductance 
values Gth- As an exception, we find peak B reproducibly 
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FIG. 3: Experimental conductance histogram data at T = 
30K fitted with lorentzian peaks. The dashed lines show the 
individual peaks. Inset: A blow-up of the peak near G = 9 Go 
fitted with two additional peaks of lower intensity is shown 
for comparison. 

at a higher conductance value than expected. This is pos- 
sibly due to size constraints in building a nanowire of two 
channels by two atoms, which for a short wire will have a 
rather large tunnelling contribution of the third channel. 

Table I lists average peak positions G exp , average peak 
widths AG, and theoretical conductance values Gth in 
columns 2 to 4. The last three columns sum up the results 
of the theoretical stability calculations, taking into ac- 
count the most stable configurations only. The compar- 
ison with the experimental data shows good agreement 
in that one can ascribe a quantized conductance value to 
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TABLE I: Experimental conductance, peak widths, theoreti- 
cal conductance, and theoretical predicted shapes of the wires 
for the most important low conductance peaks. The right 
column lists the maximum temperature T max up to which the 
wires are expected to remain stable, which is related to the 
depth of the corresponding energetic minimum. Peak G is 
decomposed into three contributions as in the inset of Fig. 3. 
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each experimental peak without ambiguity. In particu- 
lar, it is possible to predict the shape of the nanowires 
seen in conductance histograms. 

Peak D, with a conductance of 4Go, is predicted to be 
barely stable in the present analysis. However, one might 
expect to find more stable nanowires with cross sections 
having a higher-order angular symmetry [40] . Such defor- 
mations cost significantly more surface energy, and can 
thus only be energetically favorable at low conductance. 
Therefore they are not likely to modify the stability dia- 
gram of Ref. [14] significantly, but one might in particular 
expect a wire with a four-fold angular symmetry to be 
stable at a conductance of 4Go, which could explain the 
observation of peak D. 

The relative weight (intensity) of each peak in the ex- 
perimental histograms can also be extracted from the fit, 
and is shown in Fig. 4 at four different temperatures. 
The contribution of a stable wire to a histogram depends 
mostly on two conditions: (i) The wire has to be formed 
often enough to be statistically relevant; (ii) It needs to 
have a long enough lifetime so as to be recorded. The low- 
temperature lifetime of a nanowire is expected to increase 
with the maximum temperature up to which it is pre- 
dicted to remain linearly stable. This is confirmed by re- 
cent calculations restricted to cylindrical nanowires [41]. 
Condition (i) depends on the ability of the system to 
probe new configurations, which is determined by the 
mobility of the atoms and the density of stable geome- 
tries in configuration space. Most of the very stable wires 
being cylindrical, one might expect that highly deformed 
wires have a low probability to form. This affects mainly 
the elliptic wire (B) with an aspect ratio of e = 1.6. The 
elliptical wires (E) and (G) with a lower aspect ratio have 
a remarkably higher weight. Interestingly, peak (E) has 
a higher intensity than the neighbouring cylindrical peak 
(F). The big relative weight of peak (G) is due to an over- 
lay of contributions of an elliptical wire with G = 9 Go 
and two cylindrical wires with G — 8 Go and G = 10 Go 
as suggested from the inset of Fig. 3. The decrease in 
the intensity of the low conductance peaks with growing 
temperature might be explained by thermally induced 
breaking of the contact. 



CONCLUSIONS 

We have presented a detailed analysis of experimental 
conductance histograms for Na nanowires. Comparison 
with a theoretical stability analysis of cylindrical and el- 
liptical nanowires shows a one-to-one relation between 
observed peaks and theoretically stable geometries, and 
in particular allows a prediction of the cross-sectional 
shape of the wires. 

This research has been supported by the DFG through 
SFB 276 and the EU Network DIENOW. JB and CAS 
acknowledge support from NSF grant DMR03 12028. 



£ °- 3 



0.2 





20 K 




40 K 




50 K 




70 K 




J 



FIG. 4: Relative weights of the nine individual experimental 
conductance peaks (A to I) at four different temperatures, 
coded via the gray scale. 
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